On the Center of Mass in General Relativity 
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Abstract. The classical notion of center of mass for an isolated system in 
general relativity is derived from the Hamiltonian formulation and represented 
by a flux integral at infinity. In contrast to mass and linear momentum which 
are well-defined for asymptotically flat manifolds, center of mass and angular 
momentum seem less well-understood, mainly because they appear as the lower 
order terms in the expansion of the data than those which determine mass and 
linear momentum. This article summarizes some of the recent developments 
concerning center of mass and its geometric interpretation using the constant 
mean curvature foliation near infinity. Several equivalent notions of center of 
mass are also discussed. 



1. Introduction 

General relativity is the theory of the Einstein equation. Let (M , g) be a four 
dimensional manifold with a Lorentz metric g. The Einstein equation says 

G := Ric(g) - ±R(g)g = T, 

where Ric(g) and R(g) are the Ricci tensor and the scalar curvature of g respec- 
tively. The Einstein tensor G corresponds to the gravitational field of spacetime, 
while the stress-energy-momentum tensor T is related to the matter content of 
spacetime. 

If M is a spacelike hypersurface in (M.,g), we denote the induced Riemannian 
metric by g and the induced second fundamental form by K. It follows from the 
Gauss and Codazzi equations combined with the Einstein equation that g and K 
must satisfy the following constraint equations: 

R g - \K\ 2 g + {tT g Kf = 2/i, 
dw g (K -tr g Kg) = J, 

where R g is the scalar curvature of g, and /z and J are respectively the local energy 
density and the local momentum density of matter defined by T. It is more conve- 
nient to introduce the momentum tensor n := K — (tr g K)g. Then the constraint 
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equations become: 



div 9 7T = J. 



The triple (M, g, tt) satisfying the constraint equations is called an initial data set. 
An important class of initial data sets which model isolated systems is the class of 
asymptotically flat manifolds. An initial data set (M,g,w) is called asymptotically 
flat at the decay rate q greater than 1/2, if there exists a coordinate system {x} 
outside some compact set, say -Br , such that gij(x) = Sij + 0(\x\~ q ) and ~K%j{x) = 
0(\x\~ 1 ~ q ) with the corresponding decay rates on higher derivatives. If n and J 
vanish, (M, g, tt) is called vacuum. Here we allow nonzero /i and J if they satisfy 
the decay condition (j, = 0(\x\- 2 - 2q ) and J = 0{\x\^ 2q ). 

The first major investigation into general relativity as a dynamical system was 
by Arnowitt, Deser, and Misner [TJ. Later, various people further studied the 
correct formulations of total energy (the ADM mass), linear momentum, center of 
mass, and angular momentum for asymptotically flat manifolds (e.g. [101 1191 [4]) 
and their well-definedness (e.g. [3j El [9j [TT]). 

The mass and linear momentum arc defined by the flux integrals at infinity: 



The mass and linear momentum are well-defined because they are independent of 
the chosen asymptotically flat coordinates as proven by Bartnik [3] and Chrusciel 
[5]. The positive mass theorem by Schoen-Yau |20l 122] and by Witten [23] states: 
if (M,g, tt) is asymptotically flat with the dominant energy condition /i > \J\ g , 
then m is non-negative. Moreover, m = if and only if (Af, tt) is isometric to a 
hypersurface in Minkowski spacetime. 

Regge and Teitelboim |19j observed that center of mass and angular momentum 
are not generally defined for asymptotically flat manifolds unless some parity con- 
dition at infinity is imposed (see also [2l [6] for other proposed conditions). Given 
an asymptotically flat coordinate system, we denote by y odd — (f(x) — /(— x))/2 
and y cvcn = (f(x) + /(— x))/2 the odd and even parts of a function /. They are 
defined outside Br when the coordinates {x} exist. 

Definition 1.1. (M,g,n) is asymptotically flat satisfying the RT condition if 
(M, g, 7r) is asymptotically flat at a decay rate q > 1/2 and g, n are respectively 
asymptotically even/odd, i.e. 
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and similarly on higher derivatives. When we say that (M, g) is asymptotically flat 
satisfying the RT condition, we mean (M, g) satisfying the corresponding conditions 
by letting 7r = 0. 
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While all known physical examples of asymptotically flat manifolds satisfy the 
RT condition, there exist mathematical solutions to the vacuum constraint equa- 
tions which violate the RT condition, and their center of mass and angular momen- 
tum are ill-defined [4], 113) . However, these examples may not be physical in the 
sense that the evolution of the initial data is unknown. Nevertheless, under the RT 
condition, the Hamiltonian formulation of center of mass and angular momentum 
proposed by Regge-Teitelboim |19j and by Beig-0 Murchadha [4] are defined by 




where Y(p) = x x (cross product) are the rotation vector fields. 

In addition to the Hamiltonian formulation of center of mass (|1.2p . there are 
other proposed definitions such as (14.11) and (|3.4p below. A geometric interpretation 
of center of mass uses the constant mean curvature foliation in M near infinity pro- 
posed by Huisken and Yau |15j . They used the volume preserving mean curvature 
flow to construct the constant mean curvature foliation, under the condition that 
(M,g) is strongly asymptotically flat, i.e. gij = (1 + xjrr)<Jy + 0(|a;|~ 2 ). Similar 
results were proven by Ye using the inverse function theorem [24] . Metzger gen- 
eralized the results to small perturbations of the strongly asymptotically flat data 
[17j . Those authors also proved that the foliation is unique under some conditions. 
A more general uniqueness result was obtained by Qing and Tian [18] . For strongly 
asymptotically flat metric which is conformally flat near infinity, Corvino and Wu 
proved the geometric center of Huisken- Yau's foliation is equal to center of mass 
[9] , and we later removed the condition of being conformally flat . 

However, there are various interesting physical solutions to the constraint equa- 
tions which are not strongly asymptotically flat. For example, the family of Kerr 
solutions is a family of the exact solutions to the vacuum Einstein equation which 
model the rotating black holes with angular momentum. They are not strongly 
asymptotically flat, but they do satisfy the RT condition. 

In [12] . we generalized the earlier results of the constant mean curvature fo- 
liation to asymptotically flat manifolds with the RT condition, which is a natural 
condition to impose when center of mass is discussed. Furthermore, the foliation 
that we constructed is unique under some conditions. From our construction, we 
can show that the foliation is asymptotically concentric, whose geometric center is 
consistent with the classical notion of center of mass. 

A new ingredient in constructing the foliation is that for any asymptotically flat 
manifold with the RT condition, one can find a canonical family of the approximate 
spheres S(p,R). They are constructed from perturbing the coordinate spheres 
Sr{p) ■— {x : \x — p\ = R} centered at p of radius R. The approximate spheres are 
better adapted to the asymptotic symmetry than the coordinate spheres. Although 
their mean curvature is not exactly constant, they share many nice properties with 
the constant mean curvature surfaces constructed from them. For example, the 
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leading order terms of the lowest two eigenvalues of the stability operator are —2/R 2 
and 6m/ R 3 . Moreover, {S(p, R)} form a foliation centered at a fixed p for R large 
when m > 0. In many applications, the approximate spheres should be as good as 
the constant mean curvature surfaces, except that they are insensitive to center of 
mass. 

In order to construct the constant mean curvature surfaces from the approxi- 
mate spheres, we have to choose the center p of S(p, R), and a suitable choice of p 
would allow us to find a nearby constant mean curvature surface. The dependence 
of the choice of p on center of mass is based on the new observation of the following 
identity 



where H$ is the mean curvature of Sr(p). This identity can be thought as an 
alternative definition of center of mass. It also enables us to show that the geometric 
center of each constant mean curvature surface is close to center of mass, and the 
limit of the geometric centers is precisely center of mass. 

The analytic underpinnings of the above identity rely crucially on the density 
theorem (Theorem 12.31 below) that we obtained in [11] . The density theorem is 
a refinement of the density theorem established by Corvino and Schoen [§]. An 
initial data set (M, g, n) is said to have harmonic asymptotics if, outside a compact 
set, g = u A 5 and ir = u 2 (Lx5 — (divX)S) for some function u and vector field X . 
Because of the constraint equations, u and Xi are harmonic up to the lower order 
terms. Moreover, the quantities m, P, C, J can be read off of the expansion of u 
and X. The density theorems involve finding suitable spaces so that initial data 
sets with harmonic asymptotics form a dense subset of asymptotically flat data sets 
in a weighted Sobolev topology, and the physical quantities are continuous in that 
topology. 

The method to construct the approximate sequence of the initial data sets in 
the density theorems can be employed to generate more initial data sets, often with 
specified asymptotics and physical properties. Several further applications can be 
found in, for example |13|, 114] . 

The article is organized as follows. We present the density theorems in Section 
2, in which a slightly more general proposition is proven. In Section 3, we sketch the 
construction of the constant mean curvature foliation. Several equivalent notions 
of center of mass are then discussed in Section 4. 



In the previous section, the definition of asymptotically flat manifolds involves 
the pointwise regularity on g and w. A slightly weaker regularity condition on 
asymptotically flat manifolds is defined by weighted Sobolev spaces (see more de- 
tailed discussions in, for example [3]). 

Definition 2.1 (Weighted Sobolev spaces). For a non-negative integer k, a non- 
negative real number p, and a real number q, we say / S W'l'^(M) if 




(a; 1 - p l )H s da = 8irm(p l - C l ) + 0{R 



1_a «) for 1= 1,2,3, 



2. The density theorems 
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where a is a multi-index and p is a continuous function with p 
coordinates {x} are defined. When p = oo, 

ll/ll^ r( M)= E es S sup|D<7|p W +«. 

|a|<fc 

Using the weighted Sobolev norms, we say that (M, g, tt) is asymptotically fiat 
at the decay rate q if (g - 5, tt) G W^(M) x W^f_ 9 (M), and in addition (M, 5, tt) 
satisfies the RT condition if ( 5 odd , 7r ovcn ) G Wlf_ q (M \ B Ro ) x W*f_ g (M \ B Ro ). 

The proof of the positive mass theorem by Schoen and Yau |20j was originally 
stated for strongly asymptotically flat data. Later, they extended the previous 
proof to allow general asymptotically flat data by a density argument |21j . They 
observed that, given a scalar-flat metric g, there exists a sequence of scalar-flat 
metrics with harmonic asymptotics whose mass converges to the mass of g. To 
generalize their density result, one has to consider not only the scalar curvature 
equation but the full constraint equations, which is more subtle. The following 
theorem by Corvino and Schoen is the analogue for the full constraint equations. 

Theorem 2.2 (Corvino-Schoen [8]). Suppose p > 3/2 and q G (1/2,1). Let 

(Sij ~ fiiji^ij) G Wjg(M) x Wjf_ q (M) be a vacuum initial data set. There is a 
sequence of solutions (gk,^k) with harmonic asymptotics. Given any e > 0, there 
exist ko > so that 

11.9 - 9k\\ W 2 j^{M) ^ e > IK - tk Wwif (M) - £ ' for a11 k - 

Moreover, the mass and linear momentum of (g~k, 7?fc) are within e of those of {g, tt). 

The theorem says that the solutions with harmonic asymptotics are dense 
among general asymptotically flat solutions. Moreover, the mass and linear mo- 
mentum, which can be explicitly expressed in the expansion of the solutions with 
harmonic asymptotics, converge to those of the original initial data set in these 
weighted Sobolev spaces. However, in the above theorem the center of mass and an- 
gular momentum may not converge, neither are they defined generally for asymptot- 
ically flat manifolds. In the following, we show that in some more refined weighted 
Sobolev spaces, solutions with harmonic asymptotics form a dense subset inside 
asymptotically flat solutions with the RT condition so that the center of mass and 
angular momentum are continuous in that topology. Also, we can remove the con- 
dition that {g, 7r) is vacuum. 

Theorem 2.3 f |llj ). Suppose p > 3/2 and q G (1/2,1). Let (gij — Sij,%ij) G 
Wl*(M)xW*'i- q (M) be an initial data set. Suppose that ( 5 ° dd , 7if/ cn ) G Wlf_ q (M\ 

B^xWlt^MXB,,,). 

There is a sequence of data (gk,^k) with harmonic asymptotics. Given any 
e > 0, (gk,TTk) is within an e-neighborhood of (g, it) in W_' q (M) x W 1 _'[_ q (M) for 
k large as in the above theorem. Moreover, there exist R and k = k (R) so that 

\\9k dd \\w 2 _f_ q (M\B R ) - C ' ll^ Ven |lwi' 2 %(M\BB) - C > fOT a11 k - k °- 

Furthermore, mass, linear momentum, center of mass, angular momentum of (gk,^k) 
are within e of those of (g, it). 

The existence statement of the sequence of solutions with harmonic asymp- 
totics in the above two theorems involves solving an elliptic system for Uk and 



= |x| when the 
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Xk by the inverse function theorem. In the proof of Theorem I2.3[ we study more 
carefully the elliptic system that Uk, Xk and u^ dd , ^£ dd satisfy. Then we apply the 
boost trap argument to improve the decay rates of w£ dd and X£ dd and prove that 
{(i?fc, Sfc)} also satisfy the RT condition. The full proofs of the above two theorems 
are rather technical. Some arguments there can prove the following proposition 
which is slightly more general (in the sense that the sequence {(gk^k)} need not 
have harmonic asymptotics) and is of independent interest. 

Proposition 2.4. Suppose that p > 1 and q > 1/2. Assume that [g, it) and (jjk, 7ffc) 
are asymptotically flat. If for any e > 0, there exits fco so that 

IIS - 9k\\ w i.p {M \R ) < e > Ik ~ ^llw^f_ 5 (M\R ) ^ e ' for a11 k ^ fc o, 

then the mass and linear momentum of (g~k,Ttk) are within e of those of (g,Tr). 

In addition to the above conditions, if (g,ir) and {(<7fc,7Tfc)} satisfy the RT con- 
dition and if there is a constant c so that for k large, 

(2-1) ll(5-3 fe ) odd || w!f _ ?(M \ S „ ) <c, \\(«-^r en \\wi?_ q( M\B Ro) <C 

then the mass, linear momentum, center of mass, angular momentum of (gk,^k) 
are within e of those of (g, 7r). 

Proof. By the definition of mass and the divergence theorem, for any s large, 
16Trm(g k ) = lim / y.[(9k)ij,i ~ (p*)«,i]rT da 



l ,3 

X j 

[(gk)ij,% - (gk)u,j] —\ da ° 



.'• -s ■ ■ 



^[{9k)ij,ij - {9k)ii,jj]dx. 



\*\>- itJ 

In the volume integral [(9k)ij,ij ~ {9k)n,jj] 1S the leading order term of the 
scalar curvature R(g~k)- Hence, by the constraint equations, the integrand is of the 
order 0(|a;| -2-29 ), so there is a constant so such that 

/ J2[(3k)ij,ij ~ (9k)ii,jj] dx < ~, for all s > s - 
J\x\>s -Q 4 

To handle the boundary term, we symbolically denote 

_ _ x j 

y , [9ij,i - 9ii,j ~ {{9k)ij,i ~ (9k)ii,j)] -T7 da 



< I \D(g - g k )\da . 

\x\=s 



First, if p > 3, by the Sobolev embedding inequality [3j (1.9)], 

sup \D(g - g k )\\x\ 1+ « < c\\D(g - g k )\\ w x,, (M , B } 

M\B Ro X R °> 

< 49 -9k\\ w 2,p {M \ BRo y 

Fix e and sq. There exists k' so that 

1 1 .9 - 9k 1 1 w 2 j?(m\b Ro ) ^ 8( J s i-i ' for a11 k ~ k '° ■ 
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Therefore, 



/ \D(g - g k )\da < max (\D(g - g k )\s 1 +q )4ns 1 - q < ~. 

J\x\=s |x|=a ^ 



Combining the above estimates, we prove that \m(g) — m(g k )\ < e for all /c > k' . 
When 1 < p < 3, by the definition of the weighted Sobolev norm, 



!°° I (\D(9 - 9 k )\r 1+ Tr- 3 da r 2 dr < \\g - ft||; 2 

J B J\x\=r -i y x °' 

Therefore, for each fc, there exists s k with s k — > oo so that 

/ (|^( 5 -^)| s i+Y^ 1 da <^ 1/2 , 

J|:r|=s.. z 



< OO. 



F|=Sfe 

for otherwise the volume integral would diverge. Then by the Holder inequality, for 
-ir = 1 — i 



M=Sfc 



l-D(.9 - 9k)\ d(T < 



< 


u 










< 





QDig-g^sl+rsfdao) a\ 



x\=s k 

2p T p* 



Notice that 



t^ + ^-— 9<0forl<p<3 and g > 1/2. Then there exists k' Q so 
that the above surface integral is less than e/2 for all k > k' . Combining the above 
estimates, we conclude \m(g) — m(g k )\ < e. The proof of the convergence of linear 
momentum is similar. 

If we have additionally condition (|2.ip . we show that center of mass converges. 
By the definition of center of mass and the divergence theorem, 



dan 



odd 



16wm(g)C l (g) 

I X ^ ] (9ij,ij ~ dii-jj) dx 

J\x\>s i} 

T((^ d ),-(^ d ),)o"Efe dd ^ 

i,j i ^ 

•/ I X I > 3 „■ „' 



The volume integral is of the order 0(s 1 ~ 2q ) by the constraint equations and con- 
dition (|2.ip . so it is less than e/4 for all s large enough. For the boundary term, we 
need to estimate the surface integral 

[\x\\D(g k -gy dd \ + \(g k - g r dd \}da . 



rfdO 
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We then proceed as above and show that, given e > 0, there is a sequence {sk} 
with sk — > oo so that for k large 

/ [\x\\D(g k ~g)° dd \ + \(gk-g)° dd \] da < 

J\x\=s,. 1 



\x\=s k 



We can then conclude that \C (g) — G (g k )\ < £• The proof for the angular momen- 
tum is analogous. □ 

The second part of the statement in the above proposition is optimal in the sense 
that condition (|2.ip is necessary. In the recent joint work with Rick Schoen and Mu 
Tao Wang [14] . we are able to construct perturbed data (<?&, fffc) close to some given 
asymptotically fiat data with the RT condition (g, %) in wl'%(M) x W^f_ q (M). The 
energy-momentum vector of the new data sets is equal to that of (g, tt), but their 
center of mass and angular momentum can be arbitrarily prescribed. In particular, 
the perturbed data sets violate condition (12.11) . 

3. The constant mean curvature foliation near infinity 

In the asymptotically flat region, the mean curvature Hs of the coordinate 
sphere Sn(p) is close to the constant 2/R, which is the mean curvature of Sn(p) 
in Euclidean space. The difference Hg — 2/R is measured by h = g — 5. Because 
Hs(x) = divg(v g ) where div g is the divergence operator of (M,g) and v g is the unit 
outward normal vector field on Sr(p) with respect to g. By direct computations, 
we have 

2 1 ( xl ~ p 1 )( x3 ~ pOO^ - v k ) 



2 1 . ( 

Hs(x) = — + - h n,h{x)- 



JJ 2^ 3 ' y ' R? 



(3.D + 2EM*) (gi " p %f - E ^f^f- 

1 ■ / s x 3 — p 3 v-^ hu (x) _ . . 



2 ^ , R ^ R 

h3 » 

where E (x) = C^iT 1 - 2 "?) and E$ dd (x) = 0{R- 2 - 2q ). 

We first examine (|3.ip when g is strongly asymptotically flat, i.e. g = (1 + 
^)<5 + 0(|x|- 2 ). Substituting h tj = ^8 VJ + 0{\x\- 2 ) into d37TJ and simplifying 
the summation terms, we derive 

„ , 2 4m &m(x — p) ■ p 9m 2 . _, . 

The linearized mean curvature operator on Sr{p) has an approximate kernel spanned 
by {a; 1 — p 1 ,^ 2 — p 2 ,x 3 — p 3 }. It was then observed by Ye [24 that the projection 
of the last three terms of the right hand side above into the approximate kernel can 
be annihilated by choosing suitable p, which we later identified with center of mass 
[11] . Then the inverse function theorem can be applied to find the constant mean 
curvature surface near the coordinate sphere Sji(p). The two surfaces are close 
in the sense that is the graph of 4> over Sr(j>) and ||0||c 2 ' a (S H (p)) i s bounded by 
a constant uniformly in R. 

In our situation, because g has weaker asymptotics and decay rates, the inverse 
function theorem cannot be directly applied because the asymptotics of H$ — 2/R 
in (|3.1|) may be far from being a constant. It simply reflects the fact that the 
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coordinate spheres are defined with respect to some asymptotically flat coordinate 
chart which is not canonical, so the constant mean curvature surfaces may not be 
close to the coordinate spheres. One may consider to choose a better asymptotically 
flat coordinate chart, but it seems hard to handle and simplify all the summation 
terms in (|3.ip . Instead, we explicitly construct a family of approximate spheres 
S (p, R) which reflect the asymptotics better than the coordinate spheres for each p 
and large R. Moreover, from our construction, the approximate spheres also adapt 
better the RT condition [12] . 

Lemma 3.1. Let (M,g) be asymptotically flat satisfying the RT condition. There 
exists a constant c independent of R so that, for each p and for R large, there is an 
approximate sphere 

S(p, R) = {x + <j>{x)v g : 4> e C 2 ' a (S R (p))} . 
Here </>* {the pull back of <fi on 5*1(0)) satisfies 

(3.2) W\\ c ^ {Sim < cR l -«, ||(0*)° dd || c2 ,. (Sl(o)) < cR-i. 
Moreover, the mean curvature of S(p, R) is 

(3.3) Hs = ^ + f + o(R-^- 2 "), 
where f := (AnR 2 )- 1 J Sr{p) f da and f = H s - 2/R. 

Remark. When q — 1, <\> is bounded by a constant, so the approximate sphere stays 
within a constant neighborhood of the coordinate sphere. However, when q < 1, 
the size of <fi may increase as R increases. Nevertheless, S(p, R) has asymptotic 
antipodal symmetry with respect to the center p as seen in (|3.2p . 

After showing the existence of the family of the approximate spheres, we prove 
that we can perturb them to construct constant mean curvature surfaces, if the 
center p is chosen correctly. The following identity (13.41) plays a key role to locate 
the center p of the approximate spheres. 

Lemma 3.2. Let (M,g) be an asymptotically flat manifold satisfying the RT con- 
dition. For I = 1,2, 3, 

(3.4) / (x l - p l ) (h s - I) do- = ^m{p l - C l ) + OiR 1 ^). 

J S R (p) V K J 

First, it is simple to see that if the metric g = u A 5, then (|3.4I) holds as follows. 
The left hand side above is 

/ ( x i- p i)( Hs -%) d* 

Js R (p) R R J 

= -p l [ V 4u 3 u j ^ dao 

Js R (p) y R 
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On the other hand, computing m and C l when g — u 4 6, we have 



Snip) 



Sr(p) 



EX J 

: 



P 



R 



da Q = -8frm + 0(R 1 - 2q ), 



£4* 



l u 3 u 



■pi 



R 



R 



dao = -SwrnC 1 + 0(R 1 - 2 ' 1 ) 



We then conclude that (|3.4j) holds for g — u 4 S. For general metric g satisfying the 
RT condition, one may tend to apply Theorem [23] at this stage. However, it seems 
that we can only prove that (|3.4j) holds for a sequence of radii {Ri} up to some 
error e. In other words, it only allows us to conclude that 



lim 

R^oo 



{x l -p L )[H, 



Sr(p) 



) dao = &Trm(p l 



C l ). 



However, for our purpose to construct the constant mean curvature surface at each 
radius R, we need to show that (|3.4[) holds for each R. Therefore, we study the 
summation terms in (|3. II) separately and observe some cancellations in the integral 
in (|3.4p after applying the divergence theorem. In particular, our argument shows 
that (|3.4I) holds when the asymptotically flat coordinates {x} are global, without 
using the density argument. When the coordinates {x} are defined only outside the 
compact set, the (outer) boundary term is then handled by the density theorem 
(Theorem 12. 3|) . The proof is computationally involved, so we omit the detailed 
computations and present only the sketch of the proof below. 



Sketch of proof. We define, for I = 1, 2, 3, 



Sr(p) 



(x l -p l ) 



{x l -p l ){x j ~p°)(x k -p k ) 



dcr - 



Because the asymptotically flat coordinates may not be defined in the interior, 
we use the divergence theorem in the annulus A = {R < \x — p\ < Using 
integration by parts and simplifying the expression, we obtain an identity containing 
purely the boundary terms 

(3.5) Z l g {Ri) ~ Z l g {R) = B l g (Ri) - B l g {R) for all R ± > R, 

where B l q {R) equals the boundary integrals: 



JS R (p) —. 



(#) " 



R 



2hij(x) 



{x l — p l ){xi — p>) 



R 3 



d<7 



Sr(p) 



-Y 

2 ^ 



hu(x)'- 



R 



+ hu(x)'- 



R 



dao 



We would like to show that T l g (R) — B l g (R) for each R large and for I = 1,2, 3. It 
suffices to prove that 



lim . l l g (Ri) 



\im B l g {Rx). 



Ri— >oo Ri— >oo 

First notice that if g = u A 5 outside a compact set, then by direct computations, for 
any R\ large (so that g = u 4 S on (p)), 

1 1 JR 1 )=B 1 JR 1 ). 
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To prove the identity for general metrics g, we can apply Theorem 12.31 There 
exists a sequence of data {(gk,^k)} with harmonic asymptotics and radii {Rk} 
with Rk — > oo. Given e > 0, there exists k' so that 

\X l g (R k )-X l §k (R k )\<e, \B l g (R k ) - B l Sk (R k )\ < e for all k > k' . 

The proof of the estimates is similar to the argument given in the proof of Propo- 
sition [2Tj Therefore, 

lim l l g (R k ) = lim B l g (R k ). 

k— >oo 3 fc— S-CXD 3 

We then conclude that I l g (R) = B l g (R) for Z = 1,2, 3. Substituting l l g (R) by B l g (R) 
into (|3.ip and ()3.4j) . and simplifying the expression, we have 

(x 1 - P 1 )(h s -^) da Q 

S R (p) V K J 



•^k(p) , ■ R 



^ / ' X 1 X 1 -p l 

) J \ — — ha — — — ) da 



0{R^ 2q ). 



is r (p) ' 

Using the definitions of mass (jl.lj) and center of mass (|1.2p . we derive p.4[) . □ 

To construct constant mean curvature surfaces from the approximate spheres 
S (p, R) , we consider the mean curvature operator 

H : C 2 ' a (S(p, R)) -> C°' a (S(p,R)) 

defined by W(ip) equal to the mean curvature of the graph of ip over S(p, R). We 
would like to construct ip £ C 2 ' a (S(p, R)) so that %{ijj) equals a constant which 
depends only on p and R. By observing (|3.3[) . it is natural to set the constant equal 
to 2/R + f. By Taylor's formula, 

H(ip) =H(Q) + A s yj + (\A\ 2 + Ric M {^ g ^ g )) 4> 

(3.6) ri 

+ / (dH(sip) -dH{0))ipds, 
Jo 

where A5 is the Laplace-Beltrami operator and A is the second fundamental form 
of S(p, R) in M, Ric M is the Ricci curvature tensor of (M 1 g), and p, g is the outward 
unit normal vector of S(p, R). Because S(p, R) is close to S P (R) in the asymptoti- 
cally flat region, 

\ A \ 2 = + 0(R- 2 - q ), \Ric M \ = 0{R- 2 - q ). 
R 1 

Therefore, by (|3 . 3[) . to solve TL(ip) = 2/R + f is equivalent to solving 

2 

(3.7) Ao(^ o $) H -(ip o $) = / — / + lower order terms, 

R z 

where Ao is the standard spherical Laplacian on Sr(j>), and the diffeomorphism 
$ : Sr(p) -)• S(p, R) is defined by <f>(x) — x + <pv g . This operator A + 2/R 2 has 
the kernel spanned by {a; 1 — p x ,x 2 —p 2 ,x 3 — p 3 }. A necessary condition to solve 
(|3.7p is that the projection of the right hand side in (|3.7j) into the kernel has to 
be zero. By Lemma 13. 2\ if m is nonzero, we can choose the center p so that the 
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projection is zero. Moreover, the center p = C + 0(R 1 ^ 2q ) is asymptotic to center 
of mass. 

Finally, (13. 7[) is then a quasi-linear elliptic equation. We can apply the Schauder 
fixed point theorem to find a solution ip and prove the following theorem. 

Theorem 3.3 ( [12] ). Assume that (M,g) is asymptotically flat satisfying the RT 
condition at the decay rate q greater than 1/2. If m is nonzero, then there exist 
surfaces {E#} with constant mean curvature Hy, r = (2 / R) + 0(R^ 1 ^ q ) for R large. 
Moreover, T,r is a CoR 1 ~ q -graph over Sr(C), i.e. 

= {x + Mx)v g : G C 2 ' a (S R (C))} 

with 

M\\c^( Sl (o)) ^ c o Rl ~ q > and IIW) odd |l^-( Sl (o)) < c R- q , 
where is the pull back of tpo on Si (0) . 

Moreover, the geometric center defined by 

L x dao 

(3.8) lim 7* , 

is equal to the center of mass (|1.2[) . 

Furthermore, the family of constant mean curvature surfaces constructed above 
form a foliation, if mass is positive. 

Theorem 3.4 (12 ). Under the same assumption as in Theorem \3.3\ if in addition 
m > 0, then each E^ is strictly stable and {E^} form a foliation. 

We first recall that a closed constant mean curvature surface E is stable if 

(3.9) / uL s uda g = [ [|Vu| 2 - {\A\ 2 + Ric M {v g ,v g ))u 2 } da g > 

for all function u with J^uda g — 0, where L% = —As — (\A\ 2 + Ric M {v g ,v g )) 
is the stability operator. Geometrically, a stable constant mean curvature surface 
minimizes the area among nearby surfaces which enclose the same volume. If the 
inequality in Q3.9p is strict, E is called strictly stable. 



Sketch of proof. To prove that each E# constructed in Theorem 13.31 is strictly 
stable, one would like to estimate the eigenvalue of the stability operator L^ R . Be- 
cause E/j is roughly spherical, Ls R has an approximate kernel spanned by functions 
whose mean values are almost zero. To show that L^ R indeed has no kernel, one 
has to estimate the term of (\A\ 2 — 2/R 2 ) + Ric M (v g , v g ) whose pointwise value is 
however undetermined by the condition on the asymptotics of g. We observe that 
by the Bochner-Lichnerowicz identity, we only need to derive the integral estimates. 
Hence, after some cancellations, the lowest eigenvalue Ai is then estimated by 

Ai > --^ / Rio M {u g , v g ) da g + 0(R- 2 - 2q ). 

Furthermore, we can relate the integral to mass and obtain 

/ Ric M {v g ,v g ) da g = / Ric M {v ,vo) da + 0(R~ 

87TTO „ . , „. 

= -- ^ + 0(R- 1 -«), 
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where the second inequality follows from an alternative definition of mass (see (|4.2j) 
in the next section). Therefore, when m is positive, the lowest eigenvalue Ai of 
L-£ R among functions with zero mean value is positive for R large, and then is 
strictly stable. 

Furthermore, Ls R is indeed invertible. It is easy to see that the lowest eigen- 
value Ao of Le among all test functions is negative by letting the test functions 
equal to constants. We then employ the beautiful analysis of Huisken and Yau |15j 
in our setting to show that the next eigenvalue is positive. Hence Ls R is invertible. 
Then we apply the inverse function theorem for the mean curvature operator and 
conclude that {E.r} vary smoothly in R and that locally £#/ is graph of v over 
when R' is close to R. To prove that they form a foliation, it remains to show that 
£_r and do not intersect if R ^ R'; that is, v is either strictly positive or neg- 
ative. To see this, we decompose v into the projection on the lowest eigenfunction 
and the orthogonal complement. We derive that the lowest eigenfunction is close 
to a nonzero constant by the De Giorgi-Nash-Moser theory. Then by the Schauder 
estimates, we prove that the norm of the orthogonal complement is relatively small. 
Hence, v is a nonzero constant up to the lower order terms 0(R~ q ). □ 

In order to make the geometric center of mass canonical and to define a polar 
coordinate system outside a compact set in M, one would prove that the constant 
mean curvature foliation constructed is globally unique. The inverse function the- 
orem in the proof above already gives uniqueness of the foliation locally near 
for each R. For a more global result, we have the following theorem. Below, we 
denote by r — min{|x| : x € N} the minimum radius of N. 

Theorem 3.5 ( |12j ). Assume that (M,g) is asymptotically flat with the RT con- 
dition at the decay rate q greater than 1/2 and m > 0. Then there exists o\ so that 
if N has the following properties: 

(1) N is topologically a sphere, 

(2) N has constant mean curvature H — Hs R for some R> o~\, 

(3) N is stable, 

5 — q 

(4) r > H~ a for some a satisfying — — - < a < 1, 
then N = 

When q = 1, the range of a coincides with the a priori estimates by Metzger 
[IT] . When q < 1, a has a more narrow range. The optimal result that one could 
hope for is a = 0; namely, the foliation is unique outside a fixed compact set. We 
remark that it is proven true for strongly asymptotically flat manifolds by Qing 
and Tian |18) . and it remains open for general asymptotics. 

4. Equivalence of center of mass 

The constant mean curvature foliation constructed in the previous section gives 
a geometric center of mass. By our construction, it is easy to see that the geometric 
center of mass (|3.8[) is equal to the classical Hamiltonian notion (jl.2l) because each 
constant mean curvature surface is roughly round centered at p = C + 0(R 1 ^ 2q ). 



After this paper was submitted to publish, we noticed an announcement on the uniqueness 
result by Ma [16] . 



14 



LAN-HSUAN HUANG 



In addition to the Hamiltonian and geometric notions of center of mass, another 
notion using conformal Killing vector fields proposed by R. Schoen is as follows: 

(4.1) C{ = -^lim f W'«-/ > ^<K- J = 1,2, 3 

\ x \— r i,j ' 

where Yrn = Y^i (\x\ 2 S h — 2x l x l *)-J~ : are Euclidean conformal Killing vector fields. 
If the metric g satisfies the RT condition, we can replace v g and da g in (|4. 1[) by 2L. 
and dao respectively. This intrinsic definition is motivated by a similar expression 
of the mass when Yqq is replaced by {— 2x 1 ^:}, the radial direction Euclidean 
conformal Killing vector field: 

(4.2) mj = ±- hm / £ (< " \ R ^) H**')^. 

./|x|-r 4j . V / 

Notice that when g is asymptotically flat, we can replace v g and dcr s in (|4.2[) by 2_ 
and dao respectively. The Euclidean conformal Killing vector fields {— 2x l -^~} and 

generate dilation and translation near infinity. A detailed discussion on these 
conformal Killing vector fields can be found in, for example [9 . 

One way to prove (|4.1[) and (|4.2I) is to apply the density theorems in Section 2. 
It seems that the proof of m = mi has not been published anywhere, so we present 
the proof below. 

Proposition 4.1. The notion of mass (|4.2[) is equal to the classical definition (|l.ll) . 
i.e. 

mi = m. 

Similarly, (|4.ip is equal to the classical notion of center of mass (|1.2p . i.e., 

C\ = C\ for/ = 1,2,3. 

Proof. To prove m = mi, we first consider the metric g where g = u 4 S 
outside a compact set. Because u is asymptotically harmonic, its expansion is 
1 T^T ®(\ x \~ 2 ) f° r some constant A. From (jl.lj) and direct computations, 
A = m/2. Then, the expansion of the conformal factor u is 

u = i + ^ + o(M-») 

„ -rna; 4 _ _ mSa 3mx l x j ._ 4 , 

2|x|' 3 2|x| J 2|x| D 

On the other hand, because g is conformal to the Euclidean metric, the Ricci 
curvature of g is 

Rij = -2u~ 1 d i d J u + 6u- 2 diudjU - 2(u _1 Au + u~ 2 |Vu| 2 )<%. 

Substituting the expansion of u into the above formula, we get 

jz«(-2zV = nj + 0(M" 3 )- 

\x\ 

,-U., — n/IJ-4 



Also notice that i? = -4u _1 Ait = 0(|x|" 4 ). Therefore, by (l4~2]) 

la +0(W 



to/ = — — lim 

l07T r->-oc 



Lr 



dao = m - 
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For the asymptotically flat metric g with general asymptotics, we can apply the 
the density argument by Schoen and Yau |21j (or alternatively Theorem I2.2|) : Let 
{fjk} be a sequence of metrics with harmonic asymptotics. We have 

lim m(g k ) = m(g). 

k^oo 

The same argument in the proof of Proposition 12.41 would also imply 

lim mi(g~k) = mi(g). 

k— >oo 

Because m(gk) = mi{gk) as computed above, we complete the proof. 

The proof of C l — C\ can be proceeded similarly (see Theorem 1]). □ 

To see that the expression (|4.1[) is intrinsic, we prove that the integral surfaces 
can be replaced by more general surfaces, which are defined geometrically 

Proposition 4.2. Suppose that (M,g) is asymptotically flat with the RT condition. 
Let {-Dfcjfe^i C M be closed sets such that the sets Sk = dDk are connected C 1 
surfaces without boundary which satisfy 

(1) rk — inf{|x| : x G Sk} — > oo as k — > oo, 

(2) r^ 2 area(5*fe) is bounded as k — > oo, 

(3) vo\{D k \D-} = 0(rl~), 

where = D k D {— Dj~} and 3 _ is a number less than 3. 
Then the center of mass defined by 

is independent of the sequence {Sk}- 

The first two conditions (1) and (2) on Sk are the conditions also considered 
by Bartnik 3 to ensure the mass is well-defined. The volume growth condition 
(3) allows us to consider a general class of surfaces which are roughly symmetric; 
namely, the non-symmetric region Dk \ of Dk has the volume growth slightly 
less than that of arbitrary regions in M. 

Last, we can prove that the center of mass is well-defined for a fixed time-slice. 
That is, if we compute center of mass with respect to two different asymptotically 
flat charts, say {x} and {y}, the center of mass transforms properly. Assume that F 
is the transition function between {x} and {y} and y = F(x). Bartnik proved that 
the only possible coordinate changes for asymptotically flat manifolds at infinity are 
rotation and translation 0] Corollary 3.2]. More precisely, there is a rigid motion 
of R 3 , (0},a) G 0(3, R) x M so that 

\F(x) - (Ox + a) | G Wt°°(R 3 \ B Ro ). 

Theorem 4.3. Let {x} and {y} be two asymptotically flat coordinate charts of 
(M,g) with the RT condition as above. Assume that Ci^ x and Ci tV are the centers 
of mass defined by (|4.f |) with respect to {x} and {y} respectively. Then 



Ci jV = OCi }X + a. 
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An interesting phenomenon is that compared to rotation, translation is a more 
subtle rigid motion. If the translation vector a is zero, the above formula follows 
directly from the definition of Cj and Proposition @3J If a is nonzero, Theorem l2.3l 
is involved in the proof. 
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